Detecting QGP with Charge Transfer Fluctuations 
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In this study, we analyze the recently proposed charge transfer fluctuations within a finite pseudo- 
rapidity space. As the charge transfer fluctuation is a measure of the local charge correlation 
length, it is capable of detecting inhomogeneity in the hot and dense matter created by heavy ion 
collisions. We predict that going from peripheral to central collisions, the charge transfer fluctuations 
at midrapidity should decrease substantially while the charge transfer fluctuations at the edges of 
the observation window should decrease by a small amount. These are consequences of having a 
strongly inhomogeneous matter where the QGP component is concentrated around midrapidity. We 
also show how to constrain the values of the charge correlations lengths in both the hadronic phase 
and the QGP phase using the charge transfer fluctuations. 

PACS numbers: 24.60.-k, 25.75.-q, 12.38.Mh. 



I. INTRODUCTION 

More than 20 years ago, Bjorken in his seminal paperQ 
considered the possibihty that the central plateau around 
midrapidity could be due to a hot and dense matter un- 
dergoing a boost-invariant expansion. At high enough 
collision energies, the temperature and density would be 
high enough for the created matter to be composed of 
dcconfincd quarks and giuons (quark-gluon plasma or 
QGP). 

Recent studies at RHIC have shown that extreme hot 
and dense matter has indeed been created around midra- 
pidity with the energy density well above the expected 
transition density H, |^ |M IS- However, there are 
also evidences that the boost-invariance may not be a 
feature of the created system even within the appar- 
ent plateau region. For instance, the elliptic flow mea- 
sured by PHOBOS collaboration 7] shows no discernable 
plateau around the central (pseudo-)rapidity. 

Taken together, the above evidences can be regarded 
as an indication that the spatial extent (in the pseudo- 
rapidity space) of the created QGP may be only a fraction 
of the size of the plateau region. For instance, the QGP 
component may be concentrated around the midrapid- 
ity as schematically shown in Fig^ while the rest of the 
system is mostly hadronic. 

In view of such a possibility of having an inhomoge- 
neous matter, we should ask different questions about 
the produced matter at RHIC. Namely, instead of ask- 
ing whether we have created a QGP, we should ask what 
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FIG. 1: Illustration of a possible spatial inhomogeneity of the 
matter produced in central heavy-ion reactions at RHIC. The 
full line represent all charged particles produced in the heavy- 
ion collisions. The dot-dashed line indicates the fraction of 
charged particles originating from a QGP and the dashed line 
is for the charged particles that never went through a decon- 
fined phase. At midrapidity, the QGP concentration is high 
while in the forward direction, almost all particles come from 
the conflned hadronic matter. 



fraction of the produced matter went through the decon- 
fined phase and how big was the size of the deconfined 
phase. If inhomogeneity is strong, these questions should 
be answerable by some judicious choice of observables. 

In a previous paper |^ , we have proposed a new observ- 
able, the charge transfer fluctuations, for measuring the 
local charge correlation length. Since the charge correla- 
tion lengths of a QGP and a hadronic gas can be signifi- 
cantly different ^ El, i3j , the chan ges in charge 
transfer fluctuations should signal the presence and the 
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extent of the inhomogeneity. In Ref. |a| , we considered an 
ideal case assuming nearly An detector with 100% effi- 
ciency. We argued that with such an ideal detector, the 
charge transfer fluctuations should show a local minimum 
where the QGP was formed. In this foUowup paper, we 
concentrate on a more realistic scenario. The goal is to 
predict what experiments at RHIC, STAR in particular, 
should observe. 

The charge transfer fluctuation is deflned by 0, 



(1) 



The charge transfer u{ri) is in turn deflned by the 
forward-backward charge difference: 



^iv) = [Qf{v)~Qb{v)]/2, 



(2) 



where Qpifl) is the net charge in the forward pseudo- 
rapidity (or rapidity) region of rj and QBiji) is the net 
charge in the backward pseudo-rapidity (or rapidity) re- 
gion of rj. The fluctuation Du{vi) is then a measure of the 
correlation between the charges in the forward and the 
backward regions that are separated by the cut located at 
7j. As the net charge is experimentally easier to measure 
in the pseudo-rapidity space, from now on wc will use the 
pseudo-rapidity exclusively. However, all formalism can 
be directly translated to the rapidity. 

In Ref. (3, we considered the case where the experi- 
mental pseudo-rapidity coverage is much larger than the 
extent of the QGP region. Using a single-component and 
a two-component neutral cluster models, we showed that 
the presence of a QGP component results in a local min- 
imum for Du{rj) at the location of the highest concen- 
tration of the QGP because it has a shorter charge cor- 
relation length. The size of the dip can be then used 
to infer the size of the QGP region. Data from pp and 
K^p collisions in the energy ranges of pi^h = 16GeV to 
Piab — 205 GeV show that the charge transfer fluctua- 
tions is independent of rapidity 0, lla. T his is also true 
for HIJINGlll El [13 and UrQMDlli |2l| events for 
Au-Au collisions as shown in Fig|2 

If the detector coverage is comparable or smaller than 
the QGP size, then it is not likely that the local minimum 
can be observed. Among the four experiments currently 
operating at RHIC, only the STAR detector has enough 
coverage and charge-identification capability to carry out 
the charge transfer fluctuation studies. Still, the STAR 
pseudo-rapidity coverage (\rj\ < 1) is comparable to the 
extent of the QGP region we estimated in Rcf. 8] . Hence 
a quantitative analysis is necessary to show the potential 
of the charge transfer fluctuation measurement. 

The strength of the observed signal depends critically 
on the difference between the charge correlation lengths 
in the QGP and the hadronic phases. One way to es- 
timate the difference is through the net charge fluctua- 
tions. In Refs. [2^ l23| . two teams independently showed 
that the net charge fluctuations per charged degrees of 
freedom, (AQ^) / (A^ch), can be 2 to 4 times smaller if the 
hadrons came from a hadronizing QGP rather than from 
a hot resonance gas 13]. Using neutral cluster models, 
one can show that this in turn implies that the charge 
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FIG. 2: k{rj) = (Du{r]) - {AQ^}^^ /A) /{dN,i,/dri) generated 
by UrQMD Au-Au central collisions and HIJING central col- 
lisions. The net charge fluctuation {AQ^)4Tr of the produced 
particles is in general non-zero because the net charge of the 
spectators fluctuates. Hence it must be subtracted from Du 
so that the ratio remains flnite near the beam rapidities. 



correlation length is 2 to 4 times smaller in a QGP than 
in a hadronic system. 



II. THOMAS-CHAO-QUIGG RELATIONSHIP 
AND NON-QGP MODELS 

Originally, Thomas and Quigg 13 applied the charge 
transfer fluctuations in the boost invariant case and ob- 
tained 



(3) 



where Nch/Y is the value of the boost-invariant dNch/dy 
and 7 is the charge correlation length. The proportion- 
ality constant c depends on the properties of the under- 
lying clusters. Later, Chao and Quigg generalized this 
to smooth dNch/dy cases and wrote down the following 
Thomas-Chao-Quigg relationship 15J 



Du{y) 



dNch 
dy 



(4) 



where k oc 7. They also showed that k is constant (inde- 
pendent of y) for all available elementary particle collision 
data at the time. 

These original studies used a simple neutral cluster 
model where an underlying cluster with the rapidity y 
decays into 2 charged particles and a single neutral par- 
ticle. The rapidities of these 3 decay products are given 
by (y— A, y, y+A) with A oc 7. In Ref.Q, we generalized 
this simple model so that the the joint probability den- 
sity for the charged decay products is given by (switching 
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to pseudo-rapidity) 



(5) 



where r — rj^—rj^ and Y = (i]^ +r]^)/2. Here the func- 
tion F{Y) can be interpreted as the rapidity distribution 
of the clusters and R{r) can be interpreted as the rapid- 
ity distribution of the decay products given the cluster 
rapidity Y. We then showed that the above Thomas- 
Chao-Quigg relationship Q with a constant k is exactly 
satisfied if 



i?(r) 



1 

2^ 



exp(-|r|/7) 



(6) 



with J — 2k while F{Y) is chosen so that the single 
particle distribution g(r]) = J_^dri' f{r]' ,r]) yields the 
normalized dNch/dr]. Realizing that the above R{r) is 
the Green hmction of the operator d'^/dr'^ — 1/7^, one 
obtains 



2(Afo)F(7/) 



1 



4 di]'^ 



dNa, 
drj 



(7) 



where (Mq) is the average number of the neutral clusters. 
The charge density dNch/dr] is modelled with a Wood- 
Saxon form in this study and more sophisticated fittings 
are also possible|3- 

As shown in Figl^J non-QGP models of heavy ion col- 
lisions also satisfy the Thomas-Chao-Quigg relationship 
with a constant k « 0.6—0.7 or equivalently 7 w 1.2 — 1.4. 

One should emphasize here that the Thomas-Chao- 
Quigg relationship is for the case where the whole 
(pseudo-) rapidity space is observed. If the observational 
window is limited, then the ratio 



dNch/dri 



(8) 



will no longer be independent of r/ even if 7 is constant. 
The bar over k and indicates that they are measured 
only within a finite observation window. In Ref.Q we 
showed that when the observation window is confined to 
< rjo, the charge transfer fluctuation is 



Duiv) 
where 



2(Mo) 



drj 



(AQ2) = 4(Mo) 



di]^ 



dT^+f{T^+,rj-){9) 



dr^-f{i,il) (10) 



is the net charge fluctuation within |?7| < rjo and Mq is 
the total number of the neutral clusters. 

If dNch/drj is constant within the window as is the 
case for the STAR Au-Au data at RHIC energies, it is 
clear that K{r]) still varies with rj even if k itself is con- 
stant and such variation is entirely given by the second 
term in Eq.Il. In FigOl we show HIJING, UrQMD 
and RQMD^2J, |2^ results at various centralities to- 
gether with a single component model fit. The calcu- 
lation is done with the STAR acceptance rjo = I and 
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FIG. 3; The charge transfer fluctuation is shown as a function 
of pseudo-rapidity ry. Different symbols represent standard 
centralities as shown in the symbol legend. These results are 
from 8, 400 RQMD events for ^ = 130 GeV Au-Au collisions 
and 50, 000 HIJING events at ^ = 200 GeV. STAR accep- 
tance and efliciency are taken into account in a simple way by 
restricting pT > 0.1 GeV and also randomly taking out 10% 
of the charged particles. 



Pt > 0.1 GeV. Also the STAR detection efficiency is 
taken into account in a simple way by either randomly 
taking out 10 % of charged particles (for simulations) or 
by adding 10% of uncorrelated charged particles (for 
the neutral cluster model). ^From these figures, it is 
clear that the non-QGP model results are very well de- 
scribed by a single-component neutral cluster model with 
7 « 1.75 independent of centralities. The discrepancy 
between this value of 7 and the 7 obtained in the full 
phase space study is partly due to the acceptance cuts 
and partly due to the presence of large (AQ^) within 
\v\ ^ Vo- In the present case of limited observational 
window, the net charge fluctuations should not be sub- 
tracted from Du- 



III. A QGP MODEL - CENTRAL COLLISIONS 

In RHIC energy heavy-ion reactions, the energy den- 
sity of the created systems vary with the centrality of 
the collisions. The results of HIJING simulations and its 
single-component fit should correspond to the peripheral 
collision results. In central collisions, one expects that a 
QGP is formed and is concentrated around midrapidity. 
Hence, the final state particles can have 2 different ori- 
gins in central collisions: Some particles will come from 
the hadronized QGP and others will come from the non- 
QGP hadronic component of the system. 

In Ref.jsl we used a two-component neutral cluster 
model and showed that even for such an inhomogeneous 
matter Eq.JSJ still holds if one substitutes the /(?7^,?/^) 
with the following combination of the QGP correlation 
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function /qgp and the hadronic gas correlation function 
/hg 

/(ry+, r,-) = (1 - p)/HG(ry+, r?-) + pfqcpiv^ (H) 

Here p is the fraction of the charged particles originating 
from the QGP component in the whole phase space. 

Each charge correlation function is assumed to have 
the separable form as in Eqs.lQ and In the rest 

of this manuscript, we denote the hadronic correlation 
length with 71 and the QGP correlation length with 72. 
The charge correlation length is expected to be a factor 
of 2 — 4 times smaller in the QGP than in the hadronic 
matter [T^. hence, 72 < 71- The F{Y) functions for each 
components are chosen as follows. 

(l-p)FHG(n = W{Y)-pg,{Y) (12) 
pFqgp{Y) = vg2{Y) (13) 

The functions gi and 32 are required to satisfy 

Piri) = r |^.gi((V + r/)/2)e-lV-.l/7i 
J -00 271 

^g2((V + r/)/2)e-l'''-''l/^^ (14) 
) 272 

The function W(F) is again given by the left hand side of 
Eq.Q with 7 = 71. The function p(rii) has the meaning 
of the normalized dNch/drj for the hadrons coming from 
the QGP component (the dot-dashed line in Fig^J. 
For giiY) we chose a gaussian 



9i{Y) = 



(15) 



which fixes 52(^7) to be 

7I 

52 (»7) = 



dT]' 



+ (l-^jy ^5i((r;' + r;)/2)e-l"'-"l/^^16) 

again using the fact that i?HG and i?QPG are Green func- 
tions. The resulting charged particle distribution typi- 
cally looks like FigQ] 

So far we have introduced 4 parameters, 73,71,72 and 
a\ . We can fix 



71 



1.75 



(17) 



by assuming that the HIJING data shown in Fig. O is 
consistent with peripheral collisions at RHIC^.For cen- 
tral collisions, we further fix p by requiring that it is the 
maximum possible value that satisfies -Fhg(^) > 0- This 



^ The lines shown in Fig. I^Jare consistent with 1.5 ^ 7i ^ 2.0. 
The value 71 = 1.75 provides a good fit the whole set. The 
analysis given here can be repeated with any values within this 
range with minimal changes. 



in practice means p = W{0)/gi{0). For peripheral colli- 
sions, the QGP fraction p is presumed to be zero. 

The remaining 2 parameters 72 and cti can be further 
constrained by considering the experimental data on the 
net charge fluctuations from STAR |23|. ^From Ref.|27j|. 
we can infer that for central collisions 



K{rio = 1.0) w 0.27 



(18) 



The results of numerically exploring < 72 < 1.0 
are shown in Fig. 0] At each 72, there is a single cti 
that makes k(1.0) = 0.27. With cti thus fixed, one can 
then calculate the extent of the QGP component in the 
pseudo-rapidity space by calculating p and the rms- width 



-7iV2 



(19) 



of p{ri) (c.f.Eq.lHl)). Note that ^ > 7l/^/2 = 1.24 is 
always larger than the STAR rapidity window size rjo. 

The results of the charge transfer fluctuation calcula- 
tions in the two-component model are shown Fig. ^ For 
comparison, we have also plotted results of the single- 
component model calculations in Fig. |S1 The results of 
the single-component and two-component model calcula- 
tions are distinctive enough that the measurement of R{r]) 
can clearly differentiate the two model scenarios, hence 
providing a way to show the existence of the QGP. 

For the two-component model, we find the single data 
point at 7y = 1.0 does not constrain the overall shape 
of k{t]) that much. We can find a range of possible 
parameter sets that give the same k;(1.0) — 0.27 while 
the corresponding shapes of Rij]) are all very different. 
The biggest difference among these parameter sets is the 
value of k(?7) at 77 = 0.0, where the concentration of 
the QGP component is strongest. The four lines with 
0.001 < 72 < 1.0 in Fig. 21 represent typical results in 
two-component model calculations. 

In a clear contrast, in the single component model the 
value of k(1.0) completely fixes the shape of ^(77) in the 
entire interval 0.0 < 77 < 1.0. To have k(1.0) = 0.27, the 
correlation length must be 7 = 1.3. As shown in Ref.0, 
k(?7) is proportional to 7 in the limit 7/770 ^ 1- Hence 
reduction in 7 results in the overall reduction of the K;(ry) 
in the whole range. 

^From these considerations, we can say that the mea- 
surement of the charge transfer fluctuations in the entire 
range 0.0 < 77 < 1.0 for various centralities will be a crit- 
ical test for the existence of QGP. If the cental collision 
data shows a clear reduction from the 7 = 1.3 line in 
Fig. |S1 it can only be explained by the presence of the 
second component. We also expect the amount of QGP 
would grow as one goes from peripheral to central colli- 
sions. Therefore, the most peripheral collision data for 
^(77) will behave more like the single-component results 
while the most central collision data for ^(77) will behave 
more like the two-component results. We predict that 
the most central collision data should lie between the 
two solid lines with 72 = 0.6 and 0.3 in Fig^ For com- 
parison, we have also shown two extreme cases with very 
large and very small QGP charge correlation lengthes 
where 72 = 1.0 and 0.001. 
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FIG. 4: The results of two-component model calculations 
together with the single-component result that describes the 
HIJING data. Also shown is the position of the data point 
deduced from the STAR net charge fluctuations measurement. 

In summary, we have made the following prediction for 
the charge transfer fluctuations within limited pseudo- 
rapidity range of \ri\ < 1.0. Suppose that the systems 
created by two colliding heavy ions at RHIC are inho- 
mogeneous mixtures of the QGP component and the 
non-QGP hadronic component. Further suppose that 
the QGP component is concentrated near midrapidity 
and as one goes from peripheral to central collisions, the 
amount of the created QGP increases. Then the charge 
transfer fluctuations, K{r]) = Du{r])/{dNch/dT]), should 
show the following signature of the presence of the QGP: 
Going from peripheral to central collisions, K{ri — 0.0) 
should decrease substantially from about 0.45 to 0.35 
while K,(r\ = 1.0) should decrease by a small amount from 
about 0.3 to 0.27. ^From the amount of the reductions, 
we can then infer the size and the fraction of the QGP 
matter. In contrast, if a QGP is not formed at any cen- 
tralities, then the data points from different centralities 
should all fall on the same curve. 



IV. UNCORRELATED CHARGES 

The detection efficiency for charged particles is typi- 
cally less than 100 % in real experiments. This has been 
a concern in the measurement of the net charge fluctu- 
ations. In the previous paper, we have argued that this 
effect is small and will not effect the qualitative argu- 
ments we had However, with limited observation 
window, the difference between one and two component 
model is more quantitative and the detector efficiency 
deserves some attention. 

In terms of the pair correlations, the non-ideal detec- 
tor efficiency renders some of the correlated pair to be 
uncorrelated. The relevant formulas are already worked 
out in the previous paper |^. 



FIG. 5: Single-component results. The lower four lines have 
the same k(1) values as the four two-component model results 
in Fig2] Also shown are the results for the totally uncorre- 
lated case (upper most line) and the single-component model 
result that has k;(1) = 0.27 (thin hne connecting the diamond 
symbols) . 



If the detector efficiency is e and all the charged par- 
ticles are correlated, then 1 — e of detected particles will 
become uncorrelated because their partners are not de- 
tected. The corresponding charge transfer fluctuations 
are k = eKcon- + (1 ~ e)'«uncorr, where Kcorr is for the fully 
correlated charged particles and Kuncorr is for the uncor- 
related charged particles. 

The shape of Kuncorr is shown in FigEl The shape is 
flatter than the correlated cases and it is always above 
'^coiT- The detector efficiency in STAR experiment is 
about 85 — 90%. Having 10% of uncorrelated charges 
will slightly increase K(r\) and make the overall shape a 
little bit flatter. This, however, will not change much of 
the results for k{j]). The signature for the appearance of 
the second QGP-like component is still present, and we 
can still put an upper limit on the 7qgp if a significant 
reduction of k(O.O) is observed. All the results presented 
in this paper already considered the effect of uncorrelated 
charged particles. 



V. SUMMARY 

In this paper, we proposed the charge transfer fluctu- 
ation as a good observable capable of detecting the local 
presence of a QGP in a limited pseudo-rapidity space. 
In contrast, the net charge fluctuations [2^ |23| and the 
width of the balance functionj^ are only sensitive to the 
presence of a QGP when the fraction of the QGP com- 
ponent is substantial in the whole observational window. 
Since longitudinal inhomogeneity is expected from both 
theoretical considerations and experimental observations, 
it is important to have an observable that is sensitive to 
it. Furthermore, such inhomogeneity may explain why 
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the net charge fluctuations did not show a strong signal 
even though the underlying net charge fluctuations could 
be strongly reduced in the QGP phase. 

In this study, we showed that the three hadronic mod- 
els, HIJING, RQMD and UrQMD arc consistent with 
a single-component model with fixed charge correlation 
length of about about 7 = 1.75. If a QGP is created in 
central heavy ion collisions, a significant deviation from 
this behavior is expected in the data. Specifically, if the 
QGP component is concentrated around the midrapid- 
ity and tapers off going away from midrapidity, then one 
should see the following trends in the data: 

(i) The overall values of K{r]) must decrease going from 
peripheral collisions to central collisions. This in- 
dicates that more QGP is being created. 

(ii) The value of k(1.0) should change moderately from 
around 0.30 to 0.27 from peripheral to central col- 
lisions. These values correspond to the HIJING, 
UrQMD and RQMD value for peripheral Au-Au 
collisions, and the measured value of {AQ^) / (Nch) 
for central Au-Au collision from STAR. This small 
reduction indicates that near rj = ±1, the contri- 
bution from the QGP component is already much 
reduced. 

(iii) The reduction in the value of k(O.O) should be 
larger than the reduction in the value of k(1.0) as 
the QGP component is more concentrated around 
midrapidity. Going from peripheral to central 
events, the value of k(O.O) should vary from around 
0.45 (HIJING, UrQMD, RQMD) down to 0.35. 



The value of k(O.O) puts a severe constraint on the 
value of the charge correlation length for the QGP 
component. 

The change in the value of k(O.O) may not seem large. 
But keep in mind that the value of k(1.0) cannot be lower 
than the already measured value of 0.27. Hence k(O.O), 
too. cannot go lower than that. If k(O.O) is measured to 
be close to 0.35, it is impossible to explain this without 
the presence of the second phase with a very short charge 
correlation length. 

In summary, we propose that the charge transfer fluc- 
tuation is a sensitive observable to find the presence and 
extent of the QGP created in high energy heavy ion col- 
lisions. In addition, this observable is relatively easy to 
measure and does not require the net charge conservation 
correction. We strongly suggest the experimental group 
to measure the charge transfer fluctuations. 
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